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Abstract 

The temperature behaviour of the first critical field (Bci) of superconducting thin film sam- 
ples can be determined with high accuracy using an inductive and contactless method. Driving a 
sinusoidal current in a single coil placed in front of the sample, a non zero third harmonic volt- 
age V3 is induced in it when Abrikosov vortices enter the sample. Conditions to be satisfied for 
the quantitative evaluation of Bci using this technique are detailed. As validation test, different 
type II superconductors (Nb, NbN, MgB2 and YiBa2Cu307_d under the form of thin films) have 
been measured. The comparison between experimental results, data presented in literature and 
theoretical predictions is presented and discussed. 
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I. INTRODUCTION 



The first critical field, Bci, is the field at which quanta of magnetic flux (Abrikosov 
vortices) start to penetrate in a type II superconductor, and the so called mixed state sets 
in. The most popular method to determine this field is the measurement of the first mag- 
netization cycle M(H) by using, for example, vibrating sample magnetometers (VSM) or 
superconducting quantum interference devices (SQUID)^. The effective Hqi is then deter- 
mined as the field at which M(H) departs from linearity, looking at the slope (dM / dH)=-l 
— or as the field at which the induction B is different from zenA This method is very 
effective for medium-high field values whereas it lacks in precision for superconducting ma- 
terials whose first critical field is as low as the remanent magnetic moment. In this case 
the measurement of the remanent magnetic moment itself has been proposed to determine 

An interesting and alternative way to determine the first critical field is the study of 
high order harmonics in the spectrum of the electrodynamic response of a superconductor 
to an alternating magnetic field excitation. This method, suggested years agc£, is based on 
the superposition of a small sinusoidal field on a collinear steady field and has been also 
proposed to determine the irreversibility line (Tj rr }^^, below which magnetization becomes 
irreversibile as a result of flux pinning. However, the limits of validity of this technique to 
give a quantitative and reliable measurement of B C \ has never been discussed in detail. 

As a general feature, the generation of odd harmonics is predicted by the Bean critical 
state modeR This model assumes that the supercurrent density circulating in the sample 
does not exceed a critical value that is independent from the applied magnetic field whereas 
it is a function of the Bean critical field H*, beyond which all the sample volume is pen- 
etrated by the magnetic field under the form of Abrikosov vortices. It is well known that 
this model does not take into account the existence of the lower critical field. In the critical 
state the magnetization is irreversible, which means that losses are present and proportional 
to the area enclosed by the magnetization loop&2. In particular, the magnetization is pro- 
portional to the applied magnetic field whereas it does not depend explicitly on time. This 
renders the Bean model a static description of the mixed regime that completely neglects 
the flux creep contribution to losses^. The hysteretic behavior of the magnetization in the 
critical state gives rise to non zero odd harmonics in the spectrum of the electrodynamic 
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response of superconductors exposed to an ac magnetic field^. Even harmonics can appear 
too, when a dc magnetic field is superposed to an ac magnetic fielo-^i These features 
attracted a large interest in the scientific community for theoretical analysis and measure- 
ment of high order harmonics as a function of temperature, field and frequency to study 
the static^^^^^^^^^^ - 1 ^ and dynamic properties of the flux line lattice 2 ^^. A 
useful review on the vortex regime studied by high order harmonics ac susceptibility has 
been presented in ref.— . 

The aim of this paper is to consider which conditions have to be satisfied for the use 
of the third harmonic analysis as a probe of the first critical field. As validation test, we 
will show experimental results on different type II superconductors, namely Nb, NbN, MgB 2 
and YiBa2Cus07_d (YBCO), in thin film form. The way Bci values extracted by using this 
method compare with theory and experimental results obtained by other techniques will be 
also discussed. 

II. THIRD HARMONIC ANALYSIS 

We consider the case of an ac magnetic field of amplitude h and frequency / generated by 
a far source applied to a superconducting infinite cylinder with the main axis parallel to the 
field direction (demagnetizing factor N=0). As h equals Hci(T), vortices start to penetrate 
into the superconducting sample. Under these conditions a non linear power law J oc E n in 
the current-voltage curve applies, where the exponent n — > oo in the Bean limit, n > 1 in 
the flux creep regime, n — 1 in the flux flow linear regime. As a consequence, odd harmonic 
components are produced in the spectrum of the sample response signa l 10 ' 26 , when it enters 
into a region of field and temperature in the magnetic phase diagram delimited by Bci(T) 
and the irreversibility field B irr (T). In the following, we will concentrate our attention to 
the temperature, field and frequency behavior of the third harmonic component V3. In the 
Bean limit, V 3 is predicted to be proportional tc£: 

V 3 oc h 2 ■ L (1) 

By increasing the temperature in the superconducting state at fixed h and /, V3 presents a 
maximum related to the Bean's critical field H*—^-. This maximum depends on the interplay 
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between the amplitude of the applied field h and if*, whose temperature dependence follows 
the Bean critical current density (H*(T) oc Jc(T)). At low temperature and h < if* , V3 
is determined by eq. 1, therefore it increases with temperature since Jc decreases. As the 
temperature approaches the superconducting critical temperature T c , the critical current 
density is so low that the condition h > H* is easily verified. In this regime V3 oc if*, 
thus it decreases with temperature as the critical current does. By summarizing, Vs(T) 
is strictly equal to zero in the Meissner phase, assumes values different from zero with a 
bell-shape behavior in the mixed state below the irreversibility line, and it is again equal 
to zero in the flux flow and normal state regimes. The shape of the V$(T) peak is an 
important parameter to state the quality of the sample. If the zero temperature limit of the 
critical current density is low, the third harmonic peak is extremely broadened^. Moreover, 
multiple peaks can appear when regions with different values of if* are present, as it was 
shown experimentally in ref. |28|. We will refer to the left and right onset of the third 
harmonic curve as the temperatures at which V3 starts to assume values different from zero 
for rising and lowering T respectively. To better clarify these important features, we have 
presented the third harmonic temperature behaviour respect to a generic superconducting 
phase diagram in figure 1. 

The right onset of the third harmonic peak is mostly ascribed to the irreversibility 
temperature^^. It is sensitive to the frequency-^ 1 ^ 1 ^ 1 ^ 1 ^ and the amplitude of h{t)^- 
and also to the dc magnetic field intensity^ 1 ^ 1 ^ 1 ^. Since the irreversibility line is directly 
dependent on the pinning mechanism, its frequency dependence can assume different be- 
haviours in different superconducting systems as it has been put in evidence for YBCO and 
BSCCO^. 

Despite all the interest in the nonlinear electrodynamic properties of superconductors, 
the direct relation between the left onset of V 3 (T) and H C i has never been put in evidence 
properly. In order to justify this association, a discussion on the field and frequency depen- 
dence of V3 left onset is essential. As to the former, we observe from previous works that it 
decreases with increasing amplitude of the applied ac and dc magnetic field^^^^^i, as 
expected. As to the latter, we remember that V 3 (T) shows different frequency behaviours 
depending on the absence/presence of a dc magnetic field superposed to h(t). In the absence 
of dc field, and at constant h(t) amplitude and low frequency (/ < 1kHz), the Vz(T) left 
onset is /-independent^ whereas the peak position shifts towards lower temperatures and its 
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amplitude increases with increasing frequency. All these features have been experimentally 
verified up to 5KHz on YBCO thin films. For higher frequencies, Vs(T) behaviour is pre- 
dicted to be dependent on different sources of non linearity^. A further complication arises 
in the case of granular samples, where the measured first penetration field is the Josephson 
first critical field Bcu, related to the nucleation of Josephson vortices and always lower than 
the intrinsic Bci- Bcu is expected to be frequency independent up to the Josephson plasma 
frequency^, that is usually comparable or larger than 10 GH#^. If a dc field is present, by 
increasing the frequency the left onset and the peak position shift towards higher tempera- 
tures and V3 amplitude decreases^ 1 ^'^. This has been explained by taking into account 
flux creep effects^. 

For all these reasons, we focus our attention on the V3 temperature behavior in the limit 
of low frequency (/ < 1 — 5 kHz seems to be a good compromise) and in zero dc applied 
field. Under these conditions, its measurement provides information on the left side onset 
of the third harmonic response, that can be identified as due to the first vortex penetration, 
and the associated value of the applied magnetic field can be taken as the first critical field 
at that temperature. 



III. EXPERIMENTAL 



The above discussion provides us a method to evaluate the first critical field in different 
type II superconductors under the form of thin films. To measure the third harmonic 
component of the sample response to an alternate magnetic field we used the experimental 
set up shown in figure 2a. A sinusoidal current, the reference signal of a Lock-in, is amplified 
and measured by using a high precision resistor. This signal is injected in a pancake coil 
(external diameter 4 mm) located at the end of a cryogenic insert. The sample is placed in 
front of the coil at a fixed distance [hsc]- In figure 2 b a cross section of the sample-coil 
configuration is represented. Details of the coil characteristics have been presented in ref.—. 
The output signal is filtered by using a low noise notch filter to reduce the first harmonic 
component and to avoid any saturation of the Lock-in for a better signal-to-noise ratio 
during third harmonic detection. 

In order to exactly evaluate the maximum total induction field at the sample surface, we 



used the model developed in ref. |3jJ in which a superconducting infinite slab of thickness d 
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is placed at a distance hsc from a pancake coil. This model is valid only when the sample 
dimensions are larger than the external coil diameter. Instead of solving the problem in 
terms of the external magnetic field H generated by the pancake coil, it is easier to refer to 
the vector potential A since the current injected in the coil Io lies in a plane parallel to the 
sample surface. Therefore, as in a mirror image, the induced screening supercurrents will 
be localised in the sample plane. This consideration, and the assumption of a sample with 
in-plane isotropic conductivity a, imply that only the orthoradial component of the vector 
potential A expressed in cylindrical coordinates is different from zero. The total induction 
field components B r , B z and Bg on the sample surface have been exactly calculated 3 ^. Be is 
zero, B z is completely negligible, and the only component different from zero is B r , since the 
total induction field is the sum of the applied ac magnetic field and the sample magnetization. 
The maximum amplitude of B r is: 



x = » KI (2) 

where K is the geometrical coil factor, that in our experimental configuration is set to 
2 • 10 5 mT 1 (/i 5 <7=100 /im). Our sample-coil geometry gives access only to the first critical 
field in the direction orthogonal to the sample surface. Therefore, in the following we will 
refer to the measured B^i as B^. The demagnetization factor is already taken into account 
in the extraction of the total induction field at the sample surface under the fundamental 
hypotesis of infinite slab 31 . 

To validate this technique as a powerful and simple probe of the first critical field in type 
II superconducting planar samples, we have tested a number of selected thin films having 
low, medium and high critical temperature: i) a 500 nm thick Nb film prepared by ultra- 
high vacuum cathodic arc deposition 3 - 2 - on Si wafers with Tc=9.2 K; ii) a 300 nm thick NbN 
(1111) textured film grown by DC magnetron sputtering^ 3 - on (1102) — Al 2 3 with Tc=14.5 
K; Hi) a c-axis oriented 400 nm MgB 2 epitaxial film grown on (0001)4H— SiC by HPCVD 3 ^ 
with Tc=40.5 K; iv) a 700 nm thick YBCO thin film produced by thermal evaporation 
(Theva GmbH) on (100) LaAlOs with Tc=87.9 K. The surface area of all these samples is 
10x10 mm 2 . In table I a summary of the main superconducting properties of the samples 
under test is presented. Tc values have been estimated by taking the V3 right onset. This 
is correct since Tc and Tj rr merge at the limit of small fields. For Aq and £o 5 well accepted 
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average values presented in literature have been considered. 

Measurements are performed fixing the coil current amplitude and increasing the temper- 
ature below Tc after a zero field cooling. As an example, in figure 3a a set of measurements 
of V3 versus T for the Nb sample are shown. By increasing the coil current, the applied 
magnetic field and consequently the total induction field on the sample surface increases, 
producing a broadening of the V3 response. Let us focus our attention on the curve having 
the highest intensity, obtained when the maximum induction field on the sample surface is 
about 0.8 mT: one can observe that as the temperature increases, V3 remains zero until 
the first vortex penetrates the sample (left-onset). Then, it rises, reaching a maximum, 
and then decreases with a bell shape. The left onset is T(Bq 1 ), the temperature at which 
the total induction field at the sample surface (eq. [2]) is equal to the first critical field. As 
expected, when the total induction field is decreased, the left onset shifts monotonically 
towards T C (B = 0) whereas the right onset remains quite unaffected. 

The left-onset temperature is experimentally determined by looking at both the modulus 
and phase of the third harmonic signal temperature behavior (see arrows in figure 3b). In the 
amplitude versus temperature plot, the left/right onset temperature is defined as the first 
point whose amplitude exceeds the noise level floor (threshold criterion). In the specific case 
of figure 3b, this means that the sample signal third harmonic amplitude must exceed about 
0.1 fiV . The third harmonic phase has a random behavior (see the phase versus temperature 
plot) as it should be when the sample response is merely linear (in the Meissner state or in 
the flux flow regime and normal state) and no odd harmonics are present. Therefore, the 
left/right onset is defined as the temperature beyond which the third harmonic phase starts 
following a well defined pattern. The error is estimated by taking the difference between the 
V3 onset temperature and the extrapolation of the linear portion of V3 amplitude or phase 
vs temperature curve. This error never exceeds 0.5 %. Therefore, one temperature scan 
gives one point in the versus T curve. The same experiment is repeated with different 
field amplitudes for all samples under test, in order to evaluate the first critical field for a 
fairly large temperature range. 
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IV. DISCUSSION 



In figure 4 we show the measured versus T curves for all samples under test. Exper- 
imental data have been compared to the expected theoretical values calculated by using the 
following equation^ 6 .^*^: 



Box (T) = [In («) + a («)] (3) 

where k is the Ginzburg-Landau parameter, O is the flux quantum, A is the London 



penetration depth and a (k) is a small correction defined in ref. |35|. ck(k) ~ 1/2 for large 



'35,36,38 



In figure 4a the experimental data (solid circles) and the calculated values for B^ x {t = 
T/Tc) (continuous line) in the case of the Nb film are plotted. The calculation has been 
performed by using the normalized superfluid density of the BCS standard model^. For 
k, the temperature dependence reported in ref.— has been considered. The experimental 
data are in good agreement with the theoretical predictions. Moreover, they well describe 
the tiny upward curvature observed in B^^t) near the zero field critical temperature limit. 
As explained in detail below, the presence of this upward curvature could be ascribed to a 
number of reasons other than the temperature dependence of k itself^. 

For the NbN film (figure 4b), the experimental data (solid circles) are in fair agreement 
with prediction, and measurements carried out in ref.— (open squares), however near Tc only. 
The latter data have been extracted from the temperature behavior of the upper critical field 
taken by magneto-transport measurements on high quality films and using the parameters 
collected in table 1. The continuous line represents the theoretical prediction: the expected 
dependence well agrees with the experimental data down to t=0.65. However it fails to 
justify the higher values measured at lower temperatures where B^ x {t) is overestimated. 
The reason should be ascribed to the difficulty in using the assumption of an infinite slab, 
since the induction field on the sample edges cannot be neglected in comparison with the 
expected value of the first critical field. The screening is less effective than what expected 
from the model, which results in an overestimation of the current-field conversion factor. 

for the specific case of our probing coil, the total induced magnetic 



As reported in ref. 



31 



field at the sample surface is completely negligible at the sample edges only if the ratio 
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r = rg/rc > 2.5, where r$ and tq are the sample and the external coil diameters respectively. 
By increasing the intensity of the current circulating in the probing coil, the total induction 
field spatial distribution extends over larger areas on the sample surface. When it exceeds 
sample dimensions, the total induction field at the edges can be greater than the expected 
first penetration field. This is actually what happens in the case of this film for temperature 
lower than 0.65, when the currents injected in the probing coil become greater than 35 mA. 

In the case of magnesium diboride c-axis oriented film (figure 4c) , the measured Bq x [t] 
coincides with B^it) because of our sample-coil configuration (see figure2b). We can note 
an excellent agreement between the experimental data (solid circles) and the theoretical pre- 
diction (continuous line) in the full temperature range by taking into account the two-band 
BCS model normalized superfluid density of ref. |42| in the clean limit. Our experimental 
results show that the first critical field in our sample is much higher than typical values for 
thin films^ (open circles) but lower than what reported for high quality single crystals^ 1 ^ 
(open squares). 

In the case of the YBCO film (figure 4d) , experimental data are in fair agreement with the 



45l and close to the values predicted by the 



results obtained on single crystals shown in refi 
two fluid model, again near the critical temperature only. For reduced temperatures lower 
than t=0.89, the measured first critical field is much larger than the values expected from the 
theoretical prediction. Again, as in the previous case of the NbN sample, this overestimation 
can be ascribed to the fact that the hypotesis of infnite slab fails when the current injected 
in the probing coil exceeds about 70 mA. The YBCO expected first critical field is higher 
than NbN one, that likley shifts the limit of validity of the infinite slab hypothesis to higer 
values of the probing coil current. 



V. CONCLUSIONS 



Some conclusive remarks can be made for all the investigated samples. First of all, the 
presence of an upward curvature is a common feature present in all the experimental curves. 
This could be explained by taking into account that B^iT) = B^ 2 {T)-ln(n) / '(2k 2 ) . Thus, its 
presence in the upper critical field should induce the same effect in the lower critical field. 
The upward cuvature has been indeed predicted for metals exhibiting anisotropic Fermi 
surfaces^, that is the case for quite all the systems under study. However, particularly in 
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multiband superconductors like MgB2 the temperature dependencies of Bci and Bc-2 can be 
very different and the Ginsburg-Landau parameter k depends on T. Besides that, in single- 
band s-wave superconductors, a downward curvature is expected based on both extensive 
measurements and calculations of Bci using the Eliashberg theory (see for example 1471 ) . 
Likely, pinning may play a role here. Secondly, the Bci values measured in our samples 
are in good agreement with data reported in literature for film samples of similar quality 
measured by other techniques. However, they are systematically lower than those reported 
for bulk samples. This is quite a general effect due to the reduced quality of the films in 
comparison with high quality bulk samples. 

In conclusion, we have presented a new and reliable method to measure the first critical 
field in superconducting type II materials under the form of thin film. This is done by looking 
at the third harmonic component of the response signal to an ac magnetic field excitation. 
To validate this technique we have successfully measured B^iT) in low, medium and high 
critical temperature superconductors. The main factor limiting the technique as presented 
in this paper is the film size: only samples having an area larger than the coil dimensions 
at a given coil current can be measured. In particular, when sufficiently high currents are 
injected in the probing coil, and depending on the magnitude of the expected theoretical 
first critical field, the total induction field on the sample edges should be taken into account. 
This renders the hypothesis of the sample as an infinite slab questionable. As a consequence, 
an overestimation of the current-field conversion factor and of the measured first critical field 
will result. 
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TABLE I: Summary of sample properties and parameters for the evaluation of the first critical field. 



sample 


thickness [nm] 


substrate 


Tc [K] 


A [nm] 


Co[nm] 


Nb 


500 


Si/Al 2 3 


9.2 


50 til 


40 Ui 


NbN 


300 


A1 2 3 


14.5 


^50 f33 1 41 1 48 1 49l 


5 [33.41.491 


MgB 2 


40 


SiC 


40.5 


110 \BMU 


3.5 W 


YBCO 


700 


A1 2 3 


87.9 


150 W 


1.5 i^J 



Figure captions 

Figure 1: (color online) V^(T) behaviour for an a.c. magnetic field of amplitude 
fro = h /fi for a generic type II superconductor phase diagram. V 3 (T) is strictly equal to 
zero in the Meissner phase, it assumes finite values with a bell-shape behavior in the mixed 
state below the irreversibility line, and it comes back to zero in the flux flow and normal 
state regimes. The two vertical arrows refer to the left (T(Bci)) and right (T(B irr )) onset 
respectively. 



Figure 2: a) sketch of the experimental set-up; b) details of the coil-sample configuration. 

Figure 3: (color online) procedure for the extraction of the first critical field values from 
the temperature behavior of the third harmonic component (Nb sample): a) dependence 
of the Vz left onset on the applied ac field by increasing the amplitude Iq of the current 
injected in the probing coil. The arrow indicates the increase of I ; b) details of the phase 
and modulus of the dotted curve shown in a). The arrow represents the third harmonic 
onset whereas the vertical line sets the method to determine the maximum amplitude error. 



Figure 4: (color online) Bq X measured values (solid O) as a function of the reduced 
temperature t=T/T c for the Nb (a), NbN (b), MgB 2 (c) and YBCO film (d). In (b), 



(□) data represent the experimental results reported in ref. |4lJ. In (c), (□) and (O) d ata 



represent the experimental results reported in ref. 



43 for thin films and in ref. 



36 



44 for 



sing, 
ref. 



e crystals respectively. In (d), (□) data represent the experimental results reported in 



45 



. The continuous lines in each graph describe the theoretical predictions. 
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FIGURE 2a 
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FIGURE 2b 
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FIGURE 4b 
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FIGURE 4c 
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